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Abstract 

The functional relevance of regulating proteins is often limited to specific binding sites such as 
the ends of microtubules or actin-filaments. A localization of proteins on these functional sites is 
of great importance. We present a quantitative theory for a diffusion and capture process, where 
proteins diffuse on a filament and stop diffusing when reaching the filament’s end. It is found that 
end-association after one-dimensional diffusion is the main source for tip-localization of such pro¬ 
teins. As a consequence, diffusion and capture is highly efficient in enhancing the reaction velocity 
of enzymatic reactions, where proteins and filament ends are to each other as enzyme and substrate. 

We show that the reaction velocity can effectively be described within a Michaelis-Menten frame¬ 
work. Together one-dimensional diffusion and capture beats the (three-dimensional) Smoluchowski 
diffusion limit for the rate of protein association to filament ends. 

The catalytic activity of enzymes is often restricted to specific binding sites. The ends of microtubules 
(MTs) for example are binding sites for a plethora of MT associated proteins (MAPs) (TJ. At MT ends, 
MAPs can catalyze biochemical processes J2J, or serve as substrates for other enzymes. This makes an 
efficient association of MAPs to MT tips important. Recent experiments suggest that one-dimensional 
diffusion of MAPs on MTs facilitates tip-targeting |[3]S1- This idea goes back to the concept of “reduction 
in dimensionality” suggested by Adam and Delbriick Q and has been largely applied I00- However, 
a quantitative understanding of tip-binding mediated by diffusion on the filament and subsequent capture 
at the tip has remained elusive ll3ll8 UT6l . 

Here we show that capturing at the tip is crucial for tip-localization of proteins. We present a the¬ 
ory where diffusion and capture is accurately quantified with an effective association rate constant and 
provide a result which depends only on experimentally accessible parameters. For proteins which arc en¬ 
zymatically active at filament ends, our theory predicts that diffusion and capture leads to an enhancement 
of the enzymatic reaction velocity due to stronger tip-localization. We observe that the reaction veloc¬ 
ity in dependence of the enzyme concentration closely follows a Michaelis-Menten curve and quantify 
the contribution of one-dimensional diffusion to tip-localization and enzymatic processes downstream 
thereof. 

To model the diffusive motion of proteins on a filament we consider a one-dimensional lattice of 
length l with lattice spacing a = 8.4 nm [Fig. |T}A] . The lattice corresponds to a single protofilament of a 
stabilized MT in the absence of dynamic instability. Proteins perform a random walk on the lattice with 
a hopping rate e; the diffusion constant is D = e a 2 . Each site can be occupied by only one protein; the 
system is an exclusion process ifTTll . Proteins attach to and detach from the lattice at rates cj on c and w 0 g, 
respectively, where c is the concentration of proteins in solution. The tip of the MT is represented by 
the first lattice site in our model. To account for its particular structure, different on- and off-rates are 
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assumed there, k on c and k 0 g. Proteins that bind to the tip are captured, i.e. not allowed to hop on the 
lattice, but still may detach into solution. This important condition is a critical difference between our 
model and previous approaches [ 3 l [T 6 |, see also the Supporting Material. 

A central goal of this letter is to quantify the relative contributions of diffusion and capture (tip- 
attachment after diffusion on the lattice) and end-targeting (attachment after diffusion in solution) [Fig.[lj3] 
to tip-localization. To this end we calculated the probability to find a protein at the end of a protofilament 
(the tip density p + ). In the absence of diffusion and capture, the Langmuir isotherm is obtained, 
c 


P+(c) = 


( 1 ) 


A" + c ’ 

where K = k 0 s /is the dissociation constant of the protein at the tip. However, as noted pre¬ 
viously SHI, such a model is incomplete as it does not account for the additional protein flux along 
protofilaments mediated by diffusion and capture. We have analyzed this flux by stochastic simulations 
of the model [Fig. |TJ. Surprisingly, we find that over a broad range of concentrations c, the additional 
protein current to an unoccupied reaction site J D effectively obeys first order kinetics, i.e. J D = k® n c 


[Supporting Material Fig. S2|. This observation implies that despite the complexity of the diffusion- 
reaction process one approximately retains the functional form of the Langmuir isotherm. Accounting 
for the diffusion-capture contribution to the rate of protein attachment leads to an effective dissociation 
constant 


K eS = k oS /(k on + C) • 


D ' 


( 2 ) 


We have calculated the diffusion-capture rate kff analytically, by exploiting the observed approximate 
linear reaction kinetics. We find 

u) on D/a 2 


pD = 

w off + ^Wo sD/a 2 


(3) 


Refer to the Supporting Material for a detailed derivation of Eqs. EE Together Eqs. El comprise an 
effective theory for the association of proteins to the tip which accounts for direct end-targeting as well 
as the diffusion-capture process. With Eq.[3]we are able to quantitatively predict the relative contribution 
of diffusion and capture to tip-binding for different proteins that diffuse on filaments. The results are 
shown in Fig. [2j 90 — 99% of molecules bind to the tip through one-dimensional diffusion given they 
follow diffusion and capture. 

Tip-localization due to diffusion and capture as predicted by our theory has important implications 
for enzymatically active proteins. We extended the model to investigate enzymatic reactions at the MT 
tip, where the protein-tip complex catalyzes a product at rate A: ca1 [Fig. [Ip]. In detail, we assume that 
the protein does not leave the tip after catalyzing a reaction, but only through detachment into solution. 
These model assumptions are consistent with filament polymerizing enzymes that act processively, such 
as XMAP215 for MTs |[9l[l0|, and VASP Ifl5l and formins Ifl8l for actin filaments. The assumption of a 
constant length l in our model is excellent if the rate of diffusion is fast compared to the polymerization 
rate. With the above model assumptions the reaction velocity v is determined by the tip density, v = 
p+fccat- We can apply our previous results, Eqs. Ell to obtain 

^cat^ 


v(c ) = k cat p + (c) = 


(4) 


K eS + c ' 

The above equation is reminiscent of a single-molecule Michaelis-Menten equation lfT9l l20l when X off 
is reinterpreted as Michaelis constant and substrate and enzyme concentrations are interchanged. In this 
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Figure 1: (A) Schematic of a MT with diffusive tip-binding proteins. In the bulk of 
the lattice, proteins attach to empty sites and detach. Proteins hop to neighboring 
sites but obey exclusion. At the plus-end, particles are captured. (B) Illustration of 
direct tip-attachment from solution and via diffusion and capture. (C) Proteins bind 
reversibly at the plus-end. While a protein is attached there, a reaction is catalyzed 
at rate F cat . 



0.8 


0.6 


3 g 




0.4 


O XMAP 
□ MCAK 1 
0 MCAK 2 
< CLIP-170 
O EB1 
► VASP 
V PRC1 
O NuMa 
A NCD80 
O Aurora—B 


Figure 2: The model predicts the relative contribution to tip localization of proteins 
due to diffusion and diffusion & capture (color code and solid lines), k^ n /(k^ n +k on ) 
with k on = tu on (dashed for actin: a = 6 nm). Proteins that are captured at the fil¬ 
ament end (filled symbols) and proteins where evidence for capturing is lacking 
(open symbols) are shown. Proteins that in addition have a direct enzymatic activ¬ 
ity at the filament end are XMAP215 0[1O]], MCAK 1 Q, and MCAK 2 ® on 
MTs, and VASP on actin filaments, see Ref. fT5l and personal conmiunicaton [S.D. 
Hansen and R.D. Mullins, 2014]. There are also proteins that diffuse on MTs with¬ 
out enzymatic activity at MT ends, but with roles downstream of tip-localization, 
e.g. in the protein network of MT tips (U : NCD80 ifTTll : CLIP — 170 lfl2ll : NuMa, 
PRC1, EB1 liT3l : Aurora — B Ifl4l . 
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Figure 3: (A) Comparison of the reaction velocity with (solid) and without (dashed) 
lattice diffusion and with and without capturing at the tip (circles: simulation data; 
lines: analytic results) (B) shows the reaction velocity v(c). Analytic results (lines) 
are confirmed by simulation data (circles). Parameters are L = 1000, cu 0 fj = k 0 s = 
1 sec -1 , fe ca t = 10 sec -1 cu on = k on = 0.01 sec -1 nM -1 and c = 1 nM. 


way our theory constitutes an effective Michaelis-Menten theory which accounts for end-targeting and 
diffusion and capture; instead of solving a complex many-body problem it suffices to apply a mathe¬ 
matical framework which is analogous to (single-molecule) Michaelis-Menten kinetics; the details of 
diffusion and caputure are accurately included in the effective on-rate A:®® = k on + A:/ 1 . This result is in 
accordance with experimental results for several enzymatically active proteins where Michaelis-Menten 
curves were observed for the reaction speed depending on the enzyme concentration (8] [3- Inspired by 
the processive (de)polymerase activity of (MCAK) XMAP215, we assume that enzyme and substrate are 
not decomposed in the reaction step. However, it is straightforward to include a decomposition in the 
theory: the effective dissociation constant would read K cS = (k 0 s + k ca t)/k^. 

Our analytical results, Eqs.[2]|4] agree well with simulation results of the stochastic model, as shown 
in Fig. [ 3 ] A) and (B). We find that the diffusion and capture mechanism dramatically increases A:®/ and 
thereby reduces the effective dissociation constant typically by more than one order of magnitude, e.g. 
for XMAP215 we find K eS ~ 10' ~K. In the case of long dwell times and fast diffusion e, K eS 
reduces to a particularly simple form 

K cS = (A’off/Won)/V, (5) 
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where the denominator is the square root of the average number of diffusive steps a protein performs on 
the filament. Note that one-dimensional diffusion without capturing lfl6l does not lead to a particle flux 
on the filament [Supporting Material Fig. [S4| and hence the reaction velocity is not increased [Fig.[3j\]. 
Further, the particle flux might be limited by the length of the filament: Below a threshold length l c 
(which is smaller than typical in vivo lengths of MTs) we observe a length dependent behavior of the 


reaction velocity [Supporting Material Fig. S31 where our theory is not valid. 

Our analysis reveals diffusion and capture as an efficient mechanism to circumvent the diffusion 
limit for the rate of end-targeting: Smoluchowski’s theory of three-dimensional diffusion physically 
limits the rate of direct tip-attachment from solution lf2Tl . As shown here, one-dimensional diffusion 
along a filament and subsequent capture at the filament end overcomes this limitation. This has been 
shown experimentally for MCAK Q. Our work provides an applicable theory for reaction kinetics 
facilitated by diffusion and capture: specific parameter values for diffusion, tip-association and dwell 
times can be accounted for, cf. Eqs.[3]an d0 On a broader perspective our results may also be applicable 
to other systems where one-dimensional diffusion is important @ including transcription factor binding 
on DNA 
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A Derivation of Eq. (1), (2), and (4) 


In the following we derive the reaction velocity of the model presented in the Main Text. We start by 
considering the model depicted in Fig. [ST] (A): Particles from a reservoir with concentration c can bind 
reversibly to a single reaction site with rates k on c for binding and k Q s for unbinding. As long as the 
proteins are bound to the tip, they perform a not specified reaction at rate fc cat . Inspired by the processive 
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FIGURE SI: Illustration of the chemical reaction considered in this section. Parti¬ 
cles (blue circles) from a reservoir with concentration c bind reversibly to an unoc¬ 
cupied (dashed circles) reaction site at rates k on c for binding and k D g for unbinding. 
Whilst bound, the particles catalyze a (not specified) reaction at rate fc cat . In (A) 
particles can only attach directly via the reservoir. In (B) we have an additional par¬ 
ticle flux J D to the reaction site due to a diffusion and capture mechanism described 
in the Main Text. 


(de)polymerase activity of XMAP215 (MCAK), there is no decomposition of the particle and the reaction 
site in the reaction step. Including a decomposition in the model would however be straightforward (see 
also the discussion in the Main Text). Let n + denote an occupied (n + = 1) or vacant (n+ = 0) reaction 
site. Then, the average velocity of the reaction is given by 

V = (n+)fc C at- (SI) 

Here, the average refers to an ensemble average. Since the reaction site corresponds to the last lattice site 
of one protofilament, we expect our results to be valid for experimental setups with a sufficiently large 
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and constant number of protofilaments or (due to ergodicity) for the time average with respect to a single 
protofilament. In the steady state the equation for the average occupation of the reaction site reads 

0 = 4;(n+) = k on c( 1 - (?i + )) - k oS (n +). (S2) 

at 

Solving for (n+) leads to an equation for the reaction velocity which is analogous to a single molecule 
Michaelis-Menten equation HI: 


kcri i C 

K + c 


(S3) 


with I\ = k oS /k on . Note that here c is the concentration of the enzyme, not the substrate. 

If we now have an additional particle flux J to the reaction site via one-dimensional diffusion, Fig. S1 


(B), the dynamics arc changed. Including this flux in the equation for the average reaction site occupation 
leads to 


0 = ^(ra+) = k on c( 1 - (n + )) + J - k oS {n + ) (S4) 

(1 - (n+)) - k oS (n + ). (S5) 

Note that the term J D := J/ (1—(n+)) can be interpreted as the conditional particle flux to an unoccupied 
tip. If the conditional diffusive current can be written as J D := k® n c we conserve the functional form 
of Eq. [S3] but k on is replaced by an effective on-rate for direct binding as well as diffusion and capture: 

),eff _ u I uD 
rv — rv on i ^ 0 rr 

In conclusion, observing a Michaelis-Menten curve for the reaction velocity in dependence of the 
enzyme concentration c is equivalent to the statement, that the current of particles towards an unoccupied 
reaction site due to the diffusion and capture mechanism obeys first order kinetics. In Fig. [S2| simulation 
results of our model show that the linearity condition is indeed approximately fulfilled over a broad 
parameter range. The corresponding parameter values are chosen within the typical parameter range for 
one-dimensionally diffusing proteins (see Fig. [2] in the Main Text). Note that for high concentrations 
which imply almost saturated reaction velocities we observe a, in general, non-linear behavior. In such a 
parameter range extended analytic approaches than the ones presented here would be necessary. For the 
analysis shown in this work these deviations are, however, not relevant as we are interested in quantities 
which are almost saturated at such concentrations (tip-occupation and reaction velocity; particle current 
to the tip). 


— ( k on c T 


J 


1 - (n+) 


B Derivation of Eq. (3) 

B.l Mathematical model definition 

The mathematical formulation of the model depicted in Fig. [T] of the Main Text relies on a probabilistic 
description of the lattice site occupations. A basic introduction to lattice gases and related problems can 
be found in Ref. £2]. The occupation numbers n, with i £ {0,..., L — 1} describe the configuration of 
particles on the lattice, where m = 0 and rp = 1 stand for an empty or occupied lattice site respectively. 
The equations of motion for the density on each lattice site are obtained in terms of the mean occupation 
numbers (n*) = p,. Further it is assumed that neighboring lattice sites are occupied independently, which 
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FIGURE S2: The current due to the diffusion and capture mechanism towards an 
unoccupied reaction site shows linearity in c over a broad parameter range. The 
dashed line is the concentration, where the reaction velocity is 90% saturated, v = 
0.9 k ca .f For concentrations close to the maximum reaction velocity the current 
differs from the linear behavior. This is, however, largely negligible because v and 
J are almost saturated. 


is a mean-field approximation that reads (mrii+i) = {«,){n,; + i). We will justify later on that this ap¬ 
proximation is not a restriction for the computations performed to derive k^ n . Given these preliminaries, 
the equations of motion in the bulk (sites i = 2,..., L — 1) of the lattice read Q 


dt 


Pi — ^(Pi— 1 


2 Pi T Pi-\- 1) T W on c(l Pi) i0 o SPi : 


(S6) 


where particle exclusion as well as attachment and detachment kinetics are accounted for. Particle attach¬ 
ment and detachment define the equilibrium density of particles on the lattice, also called the Langmuir 
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density, 


idnr\ C 


PLa = 


+ ^off 

Note that we consider the terminal site separately and refer to its average occupation as tip density 


(S7) 


P+ '■= Po = (no) ■ 


(S8) 


The lattice site next to the tip, i = 1, is also considered distinct from the bulk dynamics: It serves as a 
boundary for the diffusive region. The equations of motion for sites i = 0,1 are 


po = e(l - p 0 )pi + k on c(l - po) - k ofi po 


dt 

—pi = e(p 2 - pi) - epi(l - po) + w 0 nc(l - pi) - W off/ 9l 


(S9) 

(S10) 


The density in the bulk is (in leading order of a gradient expansion) governed by the following diffusion 
equation 

dt.p(x,t ) = a 2 edlp(x,t) +u; on c(l - p(x,t )) - u o{ ip(x,t ), (SI 1) 

with x = ai. This is a continuous approximation of Eq.[S6] The boundary conditions are liin J: ^ x p(x') = 
PLa and p(a) = p\. The solution is obtained in the steady state dtp(x , t) = 0 and reads 


p(x ) = p La + (pi - PLa) exp[-(x - a)/A], 


with 


A = 


D 


(SI 2) 


(S13) 


^onC T ^-'ofT 

as calculated by Klein et al. 0|. Here we used D = ea 2 . This solution is valid for x > a. At sites 
i = 0 and i = 1 the density profile is not continuous and the diffusion equation can not be applied. The 
remaining task is to determine the values for po and pi in the steady state. 


B.2 Low density approximation 


To make progress, the particle flux from site i = 2 to site i = 1 from Eq. S10 is rewritten as a derivative: 
(P 2 — pi) ~ ad x p(x)\ x=a - This corresponds again to a continuous description at the corresponding 
sites. Using Eq. |S12| we see that the (right) derivative at x = a (which is equivalent to i = 1) is 


d x p\ x =a+ = ~(Pi ~~ PLa)/A. With this result we can simplify Eq. S10 in the steady state 
Pi 


0 « — a- 


PLa 


A 


Pi(1 Po)’ 


(S14) 


where attachment and detachment rates have been assumed small. Further, note that we solve Eq. S9 
the steady state for po as a function of pi: 

Pie + k on c 


in 


P o = 


(S15) 


Pie + k 

on c + k 0 f{ 

The key relation at the basis of our theory is that tip attachment via lattice diffusion obeys first-order 
reaction kinetics. This is well confirmed by simulations over a broad and biologically relevant parameter 


region, see Fig. S2 The current to an unoccupied tip is hence approximately given by J D = A:/;' c, where 


•J) . 


k^ n is a constant, independent of the enzyme concentration. In the following, we determine the diffusive 
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current for infinitesimally low enzyme concentrations and thereby determine A 11 . In this parameter 
region correlations become negligible such that the mean-field assumption becomes valid. Further, we 
assume p\ to be small. This is well justified for low concentrations. Note also that the density along the 
lattice is minimal at i = 1. Up to first order in pi Eq.|S14|reduces to 


0 = - 


CL C. Co?onC 


sj e (w 0 ff + w on c) 


+ (o + 


ek 0 f{ 


konC T ^’ofT 


)P1- 


(S16) 


The solution of the above equation determines the tip density via Eq. S15 In our low-density approxi¬ 
mation and up to first order in c we obtain 


low-c 

P 0 


kon T (U’on c) / f^'olT T \j £ W 0 ff) 
koS 


C . 


(S17) 


B.3 Site attachment due to lattice diffusion 


The solution of p + = po allows us to determine the diffusive current. In Eq. |S17| there is an additional 
term that adds to the direct attachment rate k on which vanishes for e = 0. Identifying this term as the 
diffusive on-rate the final result reads 


k 


D 

on 


€ ^on 

Woff + V e w off 


(SI 8) 


C XMAP215 parameter values 


Model parameters 6 

s 1 

Won 

(dM s ) 1 

k^off 

s 1 

^cat 

s 1 

^off 
s 1 

kon 

(nM s) 1 

XMAP215 4.7 xlO 3 

6 xlO -5 

4.1 xlO" 1 

1.3 xlO 2 

2.6 xlO” 1 

6 xlO” 5 


TABLE SI: Parameters used in our simulation for XMAP215. The numbers were 
derived from experimental data |[5] [6]]. 


D Length dependent behavior for short filaments 

For filaments shorter than a critical length l c our model implies a length dependent tip-density, see 
Fig. [S3] Below l c the particle flux due to diffusion and capture is limited because the filament is shorter 
than the length scale of free diffusion. For XMAP215 we find l c ~ 0.5 pm. This means that for typical in 
vivo lengths of MTs, where l > l c , the enhancement of tip association through one-dimensional diffusion 
and capture is realized. However, below l c length dependent behavior might occur. 
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FIGURE S3: Below a critical length l c (dotted) the reaction velocity v depends on 
the MT length. The reaction velocity saturates at the value given by Eq. 4 in the 


Main Text. Parameters are for XMAP215, c = 300 nM, see Table SI 


E Differences to work by Helenius et al. [4] and Klein et al. 


There are several differences between our work and previous theoretical investigations of the topic. Work 
by Helenius et al. 0 did not address the question of reaction kinetics and motor occupations at the 
MT tip. Their approach based on differential equations leads to mathematical inconsistencies when 
considering a capturing mechanism. This is why we pursued a different approach using a lattice gas. It 
allows us to account for the capturing mechanism which in turn leads to an explicit expression for the tip 
density. If the tip is highly occupied, the current saturates at a value which is approximately given by the 
off-rate at the tip, k a s. This saturation effect is not included in the model by Helenius et al. 0, but is 


accounted for in our model. In Fig. S4 we compare the particle current J obtained by Helenius et al. with 
our result. For high concentrations, there are large deviations between our model and the approximation 
in 0. Work by Klein et al. 0 did not include capturing at the MT tip. Therefore there is a vanishing 


current to the tip if diffusion is fast compared to the enzymatic reaction at the tip, cf. Fig. S4 



Microtubule length l [/mi] 


FIGURE S4: Differences in the (total) protein current J to the MT tip between this 
work and previous theoretical approaches. A continuous diffusion equation with 
absorbing boundary condition leads to a linear relation J oc c 0. A lattice gas 
without capturing |0 (reflecting boundary condition) has a vanishing current to the 
tip J = 0. 
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